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ABSTRACT 


The  rate  constants  associated  with  the  deactivation  of  vibrationally 
excited  CO-,  (001)  by  collisions  with  CO  have  been  experimentally  determined 
from  300  to  900°K  by  a  laser  fluorescence  method.  These  reactions  involving 
CO  considered  were: 

kl 

COz*  (001)  +  CO  -  COz  +  CO*  (v  =  1) 
k2 

*  kco 

COz  (001)  +  CO  —  co2  +  CO 

*  k4 

CO  (v  =  1)  +  co2  -2  CO  +  co2. 

Within  experimental  error  the  rate  constant  increases  linearly 

o  3  11 

with  temperature  T  (  K)  from  a  value  of  5.  7  x  10  torr  sec  at  rocm 

temperature  to  a  value  of  11.2  x  10^  torr  *  sec  *  at  900°K.  From  500  to 

900°K  the  rate  constant  k^.^  (torr  ^  sec  ^ )  varies  with  temperature  as 

log10kCO  =  A  "  BT 

with  A  =  (>.  61  and  B  =  31.  6.  From  300  to  500°K  the  measured  values  of 
k^Q  are  greater  than  those  corresponding  to  the  cited  relation.  The  rate 
constant  k^  was  found  to  be  negligible  compared  to  k^^. 

The  probability  per  collision  of  vibrational  energy  transfer  from 
C-02  (00  1)  to  CO  (v  =  0)  was  computed  by  a  theory  involving  long  range 

forces.  The  calculated  probabilities  are  in  good  agreement  with  the 
probabilities  deduced  from  measurement  of  R  j ,  the  transfer  rate  constant 
in  the  exothermic  direction. 
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I.  INTRODUCTION 


A  C02  laser  provides  an  excellent  means  of  studying  the  collisional 
deexcitation  of  CC>2  (001).  A  pulse  of  10.6/1  radiation  from  a  Q-switched 
CO^  laser  may  be  used  to  produce  in  a  test  gas  containing  CC>2  a  nonequili¬ 
brium  population  of  CO^  (001).  ^  ^  The  return  to  equilibrium  may  be 
studied  by  monitoring  the  fluorescence  radiition  in  the  vicinity  of  4.  3/i. 

The  method  is  adaptable  to  study  in  a  temperature  range  (300  to  ~  1000°K) 
not  easily  attainable  by  shock  tube  methods  but  readily  attainable  by 
conventional  electric  heating.  A  wide  attainable  temperature  range  is 
desirable  in  order  to  establish  with  precision  the  small  temperature 
coefficients  associated  with  the  collisional  deexcitation  of  CO^  . 

In  the  absence  of  vibration- vibration  coupling,  the  fluorescense 

signal  S  decays  exponentially  with  time  following  excitation  by  a  laser 
-  Rt 

pulse,  S  =  SQe  ,  and  the  decay  constant  R,  a  function  of  temperature, 
pressure,  and  gas  composition,  may  be  determined  from  the  linear 
variation  of  log  S  with  time.  The  decay  constant  R  when  normalized  to 
unit  pressure  is  a  weighted  average  of  the  specific  rate  constants 
associated  with  the  processes 

k.  . 

*  M 

COz  (001)  +  M~  CC>2  +  M  (1) 

* 

C02  (001)  +  C02  —  CQ2  +  M  (2) 


wIutc  M  is  any  sub. stain  e  other  Ilian  CO^  present  iu  the  system.  Given  a 
knowledge  of  k(;.  '  a  value  for  l<M,  the  desired  quantity,  can  be  extracted 


I'ro.n  IT  R, 


P.  p  being  the  total  jjas  pressure, 


In  the  kinetieally  more  complicated  case  of  a  substance  M  with  a 

jjc 

vibrational  level  that  couples  to  that  of  CC>2  (001),  one  must  include  in  the 
reaction  mechanism  the  pair  of  coupling  reactions 


*  — * 
CO^'  (001)  +  M  £-  ,,  COz  +  M 

M 


as  well  as  an  additional  pair  of  degradation  reactions 


M*  f  COz  — a  M  +  COz 


M'  +  M  -*b  Mr  M.  (4b) 

In  such  cases  the  signal  S  will  vary  with  time  as  the  sum  of  two  exponentials, 
S  S[e  ^  ll  +  S2e  ^  ?},  where  both  A  and  A2  depend  on  the  rates  of  all  the 
significant  reactions.  Under  xavorable  circumstances  one  can  deduce  from 
measurement  of  A  ^ ,  A  and  the  ratio  S  /S^,  values  for  the  pertinent  rate 
constants.  Such  was  the  case  in  earlier  study  of  the  system  CC>2/N2  which 
led  to  determination  of  k^,,  k-M(N2),  k^'  (N2) ,  and  *M"(N2)  in  the  temperature 
range  from  300  to  1000°K.  The  present  study  of  the  CC>2/CO  system  was 
carried  out  to  obtain  similar  information  for  M  =  CG.  The  two  systems 
C02/N2  and  CO^/CO  differ  in  that  Eq.  (3)  for  M  =  N2  is  almost  thermo¬ 
neutral  whereas  for  M  =  CO,  the  energy  mismatch  is  about  200  cm-*.  The 
kinetic  consequences  for  k^'  and  k^' 1  of  this  mismatch  were  of  particular 
interest. 


-  C.  - 


In  a  previous  study  the  vibration-to-vibration  (V-V)  energy  transfer 

(Eq.  3)  for  CO^-N^  was  explained  by  long  range  forces  using  dipole-quadru- 
4 

pole  coupling.  The  dipole  moment  matrix  element  between  000  and  001 

states  of  CO^  interacts  with  the  quadrupole  matrix  element  between 

vibrational  states  of  having  0  and  1  quantum  of  vibrational  excitation 

leading  to  energy  transfer.  Since  then  the  theory  has  been  applied  to  a 

5-7 

number  of  situations  and  has  yielded  results  in  satisfactory  agreement 
with  experiment.  The  results  are  accurate  to  ~  25%  if  the  first  order 
theory  is  applicable  and  seems  to  be  within  a  factor  of  2  if  the  second  order 
theory  is  used.  The  present  situation  affords  an  opportunity  to  apply  the 
theory  to  a  case  where  the  energy  mismatch  is  greater  than  in  any  of  the 
cases  already  considered. 

The  experimental  values  for  the  square  of  the  matrix  element  of 

8  9 

the  dipole  moment  for  the  first  fundamental  as  well  as  for  the  permanent 

quadrupole  moment  are  available.  Calculations  for  the  variation  of  the 

quadrupole  moment10,  11  and  octupole  moment11  with  internuclear  distance 

have  been  performed.  The  molecular  parameters  needed  are  thus  available 

from  independent  sources  and  a  comparison  of  the  calculated  results  with 

the  experimental  values  of  the  V-V  transfer  rate  is  expected  to  put  the 

theory  to  a  comprehensive  test. 
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II.  EXPERIMENTAL 


A.  Apparatus 

The  experimental  apparatus,  shown  schematically  in  Fig.  1,  has 
been  described  elsewhere.  1  Consequently,  the  present  description  will  be 
less  detailed.  Laser  operating  conditions  were  chosen  to  optimize  the 
fluorescence  signals.  These  conditions  were  a  laser  pressure  of  about 
17  torr  and  a  regulated  discharge  current  of  30  ma  (12  to  13  kV)  with  a 
premixed  laser  mixture  consisting  of  11%  C02>  3%  N.,,  and  81%  He.  At 
the  usual  switching  rate  of  300  cps  a  laser  pulse  consisted  of  two  closely 
spaced  pulses  with  an  overall  length  of  about  1  1/2  fisec. 

The  reaction  cell  consisted  in  essence  of  a  1"  o.d.  stainless  steel 
tube  with  provision  for  heating  the  center  section  of  the  cell.  The  furnace, 
a  resistance  wire  type  conventional  in  design  except  for  the  adaptation  to  a 
cell  of  awkward  shape,  included  a  heavy-walled  nickel  shell  to  insure 
uniformity  of  heating.  Cell  temperature  was  measured  by  means  of  four 
24-gauge  chromel-alumel  thermocouples  in  contact  with  the  center  portion 
of  the  cell  wall.  The  cell  and  furnace  assembly  is  useable  up  to  about 

1000°K. 


The  windows  (C)  required  to  transmit  the  fluorescence  radiation 
were  of  sapphire  (3/4  inch  diameter)  brazed  with  silver  to  a  Kovar  sleeve 
(Ceramics  International  Corp.  )  which  was  itself  welded  to  the  cell  body. 
Light  emerging  from  the  cell  was  transmitted  to  the  detector  by  means  of 
a  6  inch  Vycor  light  pipe  internally  plated  with  platinum.  At  cell 
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1  ic.  I  s*  lu-m.du  diagram  «»f  flvn»r«-i»«  «  nt  r  Apparatus. 


A  '• 


temperatures  greater  than  about  700  K,  tt  wit  necessary  to  t  ap  the  pipe 
with  a  small  Irtran  2  lens  (1  unh  diameter,  I  inch  focal  length)  in  order 
to  prevent  ext  essivr  heating  of  the  detector  and  holder. 

The  detection  system  included:  (I)  a  Au:Ge  photoconductive 
detector  used  in  a  conventional  manner;  (2)  a  Keithly  I02B  amplifier 
(1000  gain,  2  cps  to  either  ISO  kc  or  1.7  Me);  (  i)  a  Tektronix  Model  S35 
oscilloscope  for  visual  monitoring  of  the  amplified  fluorescence  signal; 

(4)  a  Box-Car  Integrator  (Princeton  Applied  Research,  Model  CvV-1)  and 
an  associated  X-Y  Recorder.  An  integration  device  is  necessary  because 
of  the  very  unfavorable  ratio  of  signal  to  noise  of  single  traces  at  elevated 
cell  temperatures.  Integration  •'•ries  varied  from  S  to  20  minutes. 

The  gas-handling  system  associated  with  the  reaction  cell  included 
a  lecture  bottle  source  of  bone-dry  CO^  (J.T.  Baker,  >  99.8%  pure)  and  a 
similar  source  of  research  grade  CO  (Matheson  Co.  ).  Mixing  of  CO  and 
CO^  in  a  storage  vessel  was  accomplished  by  overnight  diffusion.  Gat, 
samples  were  admitted  as  needed  to  a  small  pyrex  manifold  with  attachments 
to  the  reaction  vessel,  to  a  Manostat  manometer  (0  -  1  atm)  and  to  a  Dubrovin 
gauge  (0  -  20  torr). 

Measurements  were  usually  carried  out  in  the  following  way.  A  gas 
sample  was  admitted  to  the  reaction  cell  and  one  or  two  X-Y  records  of  the 
box  car  signal  obtained.  Usually,  the  mixture  was  then  discarded  and 
another  sample  admitted  to  the  cell.  The  process  was  repeated  until  a 
satisfactory  set  of  records  had  been  obtained.  Periodically,  the  temperature 
of  the  four  thermocouples  in  contact  with  the  cell  wall  were  measured.  The 
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arithmetic  al  average  of  the  measured  temperatures  was  considered  to  be 
the  reaction  temperature.  The  spread  in  temperature  among  the  four  thermo 
couples  was  insignificant! 

B.  Analysis  of  Experimental  Records 

Relaxation  of  excess  vibrational  energy  in  a  mixture  of  CO^  and 
CO  involves  both  V-V  coupling  between  CO-,  and  CO  and  T-V  degradation 
with  the  consequence  that  the  observed  signal  S  is  expressible  as  the  sum 
of  two  exponential  terms. 


S  Sje'V  +  (,\L>  \2)  (5) 

The  kinetic  significance  of  the  measurable  quantities  A^,  A^,  andS^/S^ 
is  discussed  in  the  following  section.  This  section  is  devoted  primarily  to 
a  description  of  the  manner  ii  which  the  desired  data  were  obtained. 

By  proper  choice  of  the  mixture  ratio  ^AI'qq,  where  ^  and  4jcq 
are  the  mole  fractions  of  CO?  and  CO  respectively,  the  decay  constant  A 

1 

can  be  made  much  greater  than  A  In  this  study,  ratios  of  /  =  1, 

^  C'  CO 

1  2,  and  1/4  were  used.  As  a  consequence  of  the  difference  between  A^ 
and  A  the  overall  reaction  can  be  partitioned  into  two  distinct  but 
overlapping  time  periods,  an  early  time  period  during  which  Eq.  5  is 
dominated  by  the  term  S^e  ''  1*  and  a  late  time  period  dominated  by  the 
term  S^e  2^  .  The  situation  is  illustrated  in  Fig.  2. 
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TIME , t  (ARBITRARY  ORIGIN)  IN20p.sec UNITS 

Fig.  2  Schematic  of  fluorescence  record  for  CO^/CO  mixture. 
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B  1  454 


Since,  for  sufficiently  late  times  the  signal  S  is  nearly  equal  to 
,  the  decay  constant  \  ^  can  be  determined  from  the  late  time 
fluorescence  trace.  The  usual  pressure  range  used  in  obtaining  clearly 
defined  late  time  histories  was  20  to  200  torr.  The  initial  time  portion  of 
a  trace  during  which  S^e  ^  1*  contributes  to  the  signal  was  ignored  and  was 
in  most  cases  useless  because  of  distortion  resulting  from  the  inability  of 
the  measurement  system  to  respond  to  the  very  high  rate  of  change  during 
the  early  time  period.  Determination  of  A  ^  from  a  late  time  history  was 
sometimes  complicated  by  slight  base  line  shifts  (a  few  percent  of  peak 
height  in  the  direction  of  increasing  signal)  that  appeared  at  high  reaction 
temperatures  (T  >  600°K).  These  base  line  shifts  have  been  attributed  to 
a  sligh  uegree  of  gas  heating  associated  with  the  absorption  and  degradation 
of  laser  energy  and  are  discussed  elsewhere.  ^  In  such  cases  the  signal  S 
was  referred  to  the  final  base  line.  Then  A^  was  determined  in  the  usual 
manner  from  the  slope  of  the  linear  variation  of  log^S  with  time.  The 
normalized  data  A^  -  A  ^/p,  where  p  is  total  pressure,  are  shown  in 
Fig.  3  for  values  of  ^q/^qq  -  1*  l/2,  and  l/4  and  temperatures  in  the 
range  from  300  to  1000°K.  The  precision  of  measurement  decreased  from 
a  few  percent  at  room  temperature  to  about  +  10%  at  high  temperatures 
(T  >  700  K),  primarily  because  of  a  pronounced  decrease  in  the  ratio  of 
signal  to  noise  at  high  temperature. 

Once  A  has  been  determined,  additional  measurements  can  be 
made  over  the  desired  temperature  range  but  at  pressures  such  that  the 
•  arly  time  history  is  clearly  defined.  Analysis  of  the  quantity  S^e  ^  1*"  , 
a  positive  quantity  at  pressures  suitable  for  a  determination  of  A  ,  during 
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this  early  time  period  leads  to  the  decay  constant  A  ^  and  a  value  of  if 
desired.  However,  since  the  measured  signal  is  S,  the  quantity  S^e'^Z*" 
must  be  calculated  and  subtracted  from  S.  To  acc  omplish  this,  A^  =  A^p 
is  first  calculated  with  the  aid  of  Fig.  3  for  the  experimental  values  of 
temperature  and  pressure.  Then,  using  a  value  of  S,  relative  to  a  shifted 
base  line  at  high  temperature,  selected  from  the  trace  at  a  time  t  such 
that  S  ~  S^e’^Z*  ,  the  value  of  is  found.  Considering  now  the  measured 
values  of  S  in  the  early  time  period,  together  with  the  calculated  values  of 
and  A^,  the  time  dependence  of  the  function  S^e  ^  1*  is  found.  Again 
plotting  the  log  of  this  function  versus  t  leads  to  a  value  of  A  and  the 
intercept  of  the  resulting  straight  line  with  the  t  =  0  axis  determines  S^. 

In  general,  A  was  determined  over  a  range  of  5  to  20  torr  and 
the  results  averaged  and  normalized  graphically  at  low  temperatures  and 
arithmetically  at  high  temperatures.  The  normalized  data,  A  *  -  A  ^/p  are 
presented  in  Fig.  4.  Precision  of  measurement  varied  from  a  few  percent 
to  +  10%  depending  on  reaction  temperature  and  mixture  composition. 

Values  of  S^/S^  were  determined  for  =  l/4  using  a  small  absorption 

cell  filled  with  CO  placed  between  the  exit  of  the  light  pipe  and  the  detector 
in  order  to  reduce  the  contribution  of  CO  radiation  to  the  observed  signal. 
These  data,  of  limited  accuracy  because  of  the  long  extrapolation  of  the 
signal  to  the  time  origin,  are  shown  in  Table  II  along  with  numerical  results 
discussed  in  the  following  section. 

C.  Data  Analysis 

When  a  mixture  of  CO^  and  CO  is  exposed  to  a  laser  pulse, 
absorption  of  a  portion  of  the  pulse  energy  by  CO^  molecules  in  the  lower 
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NORMALIZED  DECAY  CONSTANT  X  *  (  TORR  '  SEC 


10,000 


1000 


100 


CO  /CO:  Th  t*  normalized  decay  constant  A  as  a  function  of 
temperature.  ^ 
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NORMALIZED  DECAY  CONSTANT  X.x  10' 


1.4 


Fig.  4  C02/C0:  The  normalized  decay  constant  X  j  as  a  function  of 

temperature. 
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Laser  level  (100)  creates  a  nonequilibr ium  population  of  CO^  molecules  in 
the  upper  laser  level  (00  1).  At  low  reaction  temperatures  it  may  be  assumed 
that  the  excess  vibrational  energy  is  not  redistributed  to  a  significant  degree 
throughout  the  other  levels  of  the  mode.  Similarly,  it  may  be  assumed 

that  vibrational  energy  transferred  from  CO^  (001)  to  CO  is  localized  in 

$ 

the  level  CO  (v  -  1).  For  present  purposes,  therefore,  attention  has  been 
limited  to  these  reactions: 


CO  *  (001)  +  CO  CO,  +  CO*  (v  =  1)  +  206  cm'1 
L  k2  c 


COz*  (001)  +  M  (M  =  C02,  CO)  -*3  COz  +  M 


♦  *4 

CO  (v  =  1)  +  M  (M  =  C02,  CO)  —  CO  +  M 


The  reaction  scheme  involves  four  rate  constants.  Of  these,  and  k2 
are  related  to  each  other  and  to  the  temperature  T  (°K)  by  Eq.  (9) 


K  =  e'295/T 


The  rate  constants  and  Tc^,  measures  of  the  net  rate  of  energy 
degradation,  are  weighted  averages  for  M  =  C02  and  CO  and  therefore 
vary  in  general  with  mixture  composition  as  well  as  with  temperature. 

Corresponding  to  the  reaction  mechanism  are  these  two  coupled 
linear  first  order  differential  equations  in  the  dependent  variables 
X  (CO,*  -  CC>2  )  and  y  -  (CO*  -  CO  )  where  CO  2  and  CO  are  respectively 
the  equilibrium  populations  of  CO,  (001)  and  CO*  (v  -  1). 
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•ar  =  (*kipco  -Vm*  x  +  kjPCKy 


■gf"  ■  kipco  x+  y  <’kipcK  '  ^4PM* 


It  is  known  that  both  x  and  y  vary  with  time  as 
x  =  A^-V  +  Bje-V 


a  -  A  ,  t  n  -A  -.t 
y  =  Ay  e  1  +  B?e  2 


with  decay  constants  A  and  A^  given  by 


where 


2X1,2  =  R1  ±VR1  -4R2  <X1>  X2> 


R1  =  klPCO  +  klPCK  +¥3PM  +¥4PI 


R2  '  k1^4pCOPM  +^3¥4P  M  +  k1^3PCpMK 

normalized  form  A  =  A  .  /p.  - ,  A  _  ^  =  A  -,/p. . 

1  1  M  L  2  M 


R  1  "  kl  ^CO  +  kI^CK  +  "^3  +  "^4 


R2  ”  kl"^4^CO  +"^3^4  +  k1^3^CK 


with  4<qo  anc*  eclua^  respectively  to  the  mole  fractions  of  CO  and  CO^. 
It  may  be  shown  from  the  initial  conditions  that 
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If  (liiiirtm  en»  e  radiation  (rum  (!0  doe*  not  contribute  stgntfic  antly 
to  *be  oliM'rvi'tl  nfKft.ll  S.  then  the  ratio  A  ^  H  may  be  identified  with  the 
measurable  cpiantity  S  S,.  In  that  »  ase,  F.qs.  (|7),  (|M),  and  ( 1 9)  c  tsistitute 
a  system  *>f  three  etpiatinua  in  the  three  unknowns  kj,  l  j,  and  measure- 
ment  of  t  .  *,  '  ,  * .  and  S  S,  provides  sufficient  information  to  establish 
values  for  k|,  R  and  R^.  Phis  was  the  eaac  in  an  earlier  and  aimilar 
study*  of  the  C.O,  N,  system.  Because  N  *  (v  I)  does  not  radiate,  the 
observed  signal  was  due  to  CO^*  alone  and  the  ratio  S ^ ,  S ^  was  measured 
with  sufficient  accuracy  to  demonstrate  that  k^(N^)  could  be  ignored 

T  •  * 

relative  to  k  y  In  the  present  case,  both  CO  and  CO^  can  citntribute  to 
the  signal  and  unless  radiation  from  CO  is  suppressed,  a  value  of  Sj 
is  an  unreliable  measure  of  the  ratio  Aj,  Bj.  Measurement  of  4  ^  and  \  ^ 
alone  is  then  insufficient  to  determine  kj,  Tc  ^ ,  and 

Fortunately,  the  relative  importance  of  "R  ^  and  R^  can  be  established 
using  the  t  ,  data  of  Fig.  i.  A  noteworthy  feature  of  that  data  is  that  at  a 
given  temperature,  t  is  closely  proportional  to  the  ratio  It 

i  an  be  shown  that  the  observed  proportionality  is  consistent  wtthHj  •  R^ 
but  nu  (insistent  with  the  i  oiiverse  "R  R  Consider  F.cj,  (14).  Typically, 

I 

in  this  study  R  ,  R  “  <  o.  I.  Consrcpiently,  the  square  roi*t  term  can  !>*• 

.i pproxim.ctcd  by  a  two  term  expansion  with  the  result  that 


-  I*  - 


,  I  „  -l 

1  7; 


k|VCO('  *  'ir) 


For  •  0,  Eq.  (20)  bn  omn 


k  1*  JVCK 


ki4,co( 1  * 


If,  in  addition,  "Rj  •<  k  p  as  is  probably  true,  then  Eq.  (21)  reduces  to 
Eq.  (22) 


or,  with  rearrangement  and  substitution  of  the  relation Tc^  =  k^jj^Q  +  kCO^CO’ 

2 

to  Eq.  (21) 


*  ('  *  +^r)  =  K*ckco  (' 


kc4'c 


xcrco 


In  a  similar  manner  it  can  be  shown  that  the  assumption Rj  =  0  leads  to 


Eq.  (2*) 
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Ot  tin-  two  i  iiiii|>(im'iii m  tit  k  |  ,  that  i  or responding  to  the  reaction 


CO  ♦  CO  —  CO  ¥  CO 


l  i 


is  known  to  be  negligibly  small.  Consequently ,  R^  ^4C^C  w*u>re  ^4(3 
rt'liTs  to  tin*  [irm  i‘ss 

k 


4C 

CO”  f  CO ,  —  CO  ¥  CO 


(26) 


In  terms  of  k  .  ,  Htj .  (24)  becomes 


*■  V  vco  / 


^c^c 


A 1  *  ording  to  Eq.  (2  3)  the  ()uantity  A  I  1  ¥ 


1  \  ^CO  / 


(27) 


varies  nonlinearly  with 


.....  to  a  degree  that  depends  on  the  ratio  k-Zk-^  whereas  by  Eq.  (27), 

1  /  ^C-  \ 

\  >  I  |  v  -  I  IS  strictly  proportional  to  The  different  dependence 

-  V  'co  '  c 

on  of  F.qs,  (2  3)  and  (27)  permits  a  choice  to  be  made  between  the  case 

R  ^  '^4  an(l  ^4  using  the  data  of  Fig.  3.  As  shown  in  Fig.  5  the 

1  /  VC^  \ 

variation  of  A  .  I  1  * - I  with  is  markedly  nonlinear,  indicating 

CO  j  /  \ 

t  iiat  R  .  R,.  Further,  replotting  A  (  1  4-  — -  K  as  a  function  of 

'  Mc°  / 

'V  vCO  ro8u*ts  *u  a  nearly  linear  relation,  indicating  that  "R^  is  insensitive 
to  mixture  ratio,  and  that  k^.  and  k^-Q  are  comparable  in  magnitude. 

With  the  simplification  R^  U  the  data  shown  in  Figs.  3  and  4  were 
analyzed  numeric  ally.  The  results  of  that  analysis  are  ambiguous  in  that 
t  •  a  pair  of  values  for  A  ,  there  corresponds  two  solutions  differing 
markedly  in  values  for  k  ,  R  ^  and  Aj/Bj.  The  correct  solution  can  be 

ted  either  by  physi.  al  reasoning  (R  ^  k^  at  least  at  low  temperatures) 
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BS540 


or  by  .in  approximate  im,,isumm,i  of  S^/S^.  Such  measurements  were 
available  for  ^  1,  4.  (Set*  Table  l.  )  Both  approaches  lead  to  the 

same  choice  of  solutions.  The  results  of  the  analysis  for  are  shown 
graphically  in  Fig.  6.  vVithin  experimental  error,  k^  increases  linearly 
with  temperature  unlike  the  beha\  ior  of  the  analogous  k^  for  CO^/N^ 
shown  in  Fig.  b  for  comparative  purposes. 

As  noted  earlier,  the  degradation  rate  constant  Tc^  is  a  weighted 

-j  k  T  t  .  Lr  4.  .1.  If 

c 

k. 


average  of  k^,  and  kco>  "R  ^  -  vckc  +  . 


co2  +  co2  —  co2  +  co2 


(28) 


«  CO 

co2  +  co  —  co2  +  CO  . 


(29) 


Given  a  knowledge  of  k^,,  values  of  k^.^  were  extracted  from  with  the 

results  shown  graphically  in  Fig.  7  in  conventional  Landau-Te' ler  form, 

log  | o^oo  as  a  ^uncti°n  T  ^  /t  high  temperatures  (T  >  500°K),  the 

variation  of  k^.^  with  temperature  can  be  represented  by  log^k^.Q  =  A  -  BT 

with  A  6.61  and  B  51.6.  As  in  the  case  of  k^.  ^  and  k^  ^  the  nonlinearity 

at  low  temperature  is  attributed  to  the  existence  of  a  shallow  well  in  the 

* 

potential  function  describing  the  interaction  of  C02  with  CO,  as  suggested 
by  theoretical  analysis  of  the  collisional  deactivation  of  vibrationally 


-1/3 


existed  diatomic  molecules. 


14-16 
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as  a 


TABLE  I.  Intercept  ratio  Sj/S^ 


function  of  temperature 

T(°K) 

S, 

*T» —  (Exp) 
b2 

ft) 

300 

10  -  11 

11.4 

5525 

351 

10 

9.  9 

4962 

434 

8  -  9 

8.  6 

3275 

536 

8 

7.  8 

1947 

630 

9  -  10 

7.  8 

841 

744 

7  -  8 

8.  1 

348 

858 

6  -  7 

8.  5 

213 

*  Eq.  (19),  correct  solution 
t  Eq,  (19),  wrong  solution 

=  mole  fraction  of  CO^ 
4*co  =  mole  fraction  of  CO 
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D.  I  Ih'h  rt-t  u  a  I 

The  value  i > f  tin-  specific  reaction  rate  for  vibration  to  vibration 
energy  transfer  (Fig.  (>)  eorresponds  to  a  probability  of  energy  exchange 
for  collision  P  equal  to  o.  5  x  10  ^  and  2.  I  x  10  ^  at  300  and  900°K, 
respectively.  P  increases  monotonically  with  temperature  and  at  higher 
temperatures,  where  the  signal  to  noise  problems  are  severe,  is  estimated 
to  have  a  precision  of  about  10%.  At  lower  temperatures  the  experiments 
are  more  reliable.  Interestingly  enough  the  values  of  P  obtained  by  a  Landau- 

1  —  J  g 

Teller*  or  Schwa  rt/- Slawsky- I  ler /field  type  calculation  are  of  the  magni- 

1  9 

tiule  cited  above.  For  example,  Rapp  calculates  P  for  the  reaction. 

**  .I 

A 13  (v  1)  f  CD  (v  0)  —  AB  (v  -  0)  +  CD  (v  =  1)  +  200  cm  (30) 

.  4  o 

in  the  exothermic  direction  to  be  about  8  x  10  at  300  K  and  increasing 
to  about  4  x  10  *  at  1000°K.  A  claculation  of  the  probabilities  V-V  energy 
transfer  due  to  long  range  forces  then  affords  an  opportunity  to  compare 
two  predictions  based  on  entirely  different  premises  with  each  other  and 
with  the  experimental  result.  . 

The  calculation  of  the  V-V  rates  due  to  long  range  forces  is 
divided  into  two  parts,  the  contributions  to  the  rate  in  the  first  and  second 
order  perturbation  expansions.  The  first  order  calculation  was  first 

4 

presented  by  Sharma  and  Brau  and  has  proven  its  usefulness  in  a  variety 

of  situations.  ’  1'he  second  order  calculation  to  be  used  here  has  been 

presented  1  elsewhere.  This  calculation  is  valid  when  a  .7  >  1  where  u. 

vi  vi 

is  the  energy  difference  between  the  virtual  and  initial  state  divided  by  )( 
and  7  is  the  time  duration  of  the  collision.  The  usefulness  of  this  calcula¬ 
tion  has  been  demonstrated  by  applying  it  to  the  V-V  transfer  reaction 
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C02(001)  +  N20(000)  —  C02(000)  +  N20(001)  +  125.  8  cm 


(31) 


-1 

the  rate  constant  of  which  has  been  measured^  at  room  temperature.  The 

first  order  and  second  order  calculations  are  expected  to  be  accurate  to 

about  25%  and  a  factor  of  about  2,  respectively. 

In  the  first  order  calculation  we  will  calculate  the  rate  constant 

due  to  dipole-dipole  interaction  alone.  We  will  discuss  the  contribution  of 

the  higher  multipoles  later  in  the  paper.  The  reason  for  doing  this  is  that 

only  the  dipole  moment  matrix  elements  are  known  with  confidence.  There 

are  several  calculations^  for  the  quadrupole  moment  and  all  of  these  give 

o 

reasonable  agreement  with  the  experimentally  measured  value.  Two  of 

these^’  H  calculate  quadrupole  moment  as  a  function  of  internuclear 

4 

distance.  However,  past  experience  with  N2  shows  that  a  calculation 
which  yields  a  reasonable  value  for  the  quadrupole  moment  may  not  give 
the  slope  accurately.  Since  at  this  stage  we  want  to  calculate  the  rate 
constant  using  parameters  the  value  of  which  cannot  reasonably  be  disputed, 
we  will  reserve  talking  about  the  contribution  of  higher  multipoles  until 
later. 

The  probability  of  energy  transfer  per  collision,  P,  due  to  dipole- 
dipole  coupling  is^’  ^ 

Pdd  =  p  (b  =  °)  +  i.  5  P  (b  =  d)  (32) 

where  P  (b  =  0)  and  P  (b  =  d)  are  the  values  of  P  for  impact  parameters 
equal  to  zero  and  hard  sphere  diameter,  respectively.  Previous 
experience  shows  that  P  (b  =0)  is  about  10%  of  the  total  P.  Since  our 
calculation  does  not  have  this  accuracy,  we  will  neglect  the  first  term 
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on  the  right  hand  side  of  Fq.  (52).  Following  Ref.  4  we  can  write 


dd 


-or-  <r4  n  [  ^  [-)<?>  c2(j.i  j^ooiji^x) 


■VJi 
i  -  1.2 


(33) 


where  M  is  the  reduced  mass  of  the  collision  pair,  ju ^  and  H^q  are  the 
transition  moment  for  the  v  of  CO^  (molecule  1)  and  CO  (molec  lie  2), 
respectively.  The  absolute  square  for  the  transition  dipole  moments  was 

U  3  n  jo  2 

taken  to  be  1.0  x  10  and  1.  1  x  10  (esu-cm)  for  CO^  and  CO, 

-  8 

respectively.  The  term  d  is  the  hard  sphere  diameter,  d  =  3.  82  x  10~ 

cm,^  rj  is  the  probability  of  finding  molecule  i  in  the  level  j., 

Ji  , 

C  (j^  1  j.  ;  00)  ar  the  reduced  Clebsch-Gordon  coefficients  and  (x)  is  the 

Fourier  transform  of  the  dipole -dipole  potential  given  by 


2. 


l2  (x)  -=  (0.  3333  4-  0.  2618  x")  /  (x)  4  0.  3334  x  f(|(x)  4  0.  1964  y  f  Q(x), 


x 


,  ,  /  M  \1//2  .  .  ^  n  -[  x/u  4  u^J  ,  ,  J  .  .  d  .  . 

ikT )  ’  Vx)  :£ue  u  fo (x) " '  -ar  4(x) 


(34) 


Figure  8  shows  a  plot  of  as  a  function  of  T  K.  The  rate  constant  is 
strongly  dependent  upon  temperature  and,  in  contrast  to  the  cases  previously 
investigated,  increases  with  temperature.  The  reason  for  this  is  as  follows: 
In  the  theory  of  energy  transfer  using  long  range  interactions,  certain  rota¬ 
tion  selection  rates  are  observed.  For  example,  dipole -dipole  interaction 

leads  to  A:  A  :  4  1  rotational  selection  rates.  Since  the  energy 

J  1  -2  — 

transfer  probabilities  are  large  if  u.  T  <  2,  where  u.  is  the  absolute  value 
of  the  energy  that  is  transferred  from  internal  degrees  for  freedom's 


translation  (radians  sec  and  T  is  the  time  duration  of  the  collisions, 
those  collisions  which  absorb  most  of  the  energy  mismatch  in  rotations 
are  strongly  preferred.  Since  the  energy  mismatch  that  can  be  absorbed 
into  rotations  is  proportional  to  j,  and  the  higher  rotational  levels  are 


more  populated  at  higher  temperatures,  the  cross-section  increases  with 


temperature,  the  P  v;;  T  curve  will  begin  to  level  off  at  a  temperature 

AE,  where  J _  and  J _  are  the  most 

’  m  ^  ~ 


m- 


T,  such  that  2  B  J„  +  2  J 
1  1  m  1  2  m2 

probable  values  of  the  rotation  quantum  numbers  at  the  temperature  T  and 
AE  is  the  energy  mismatch.  With  further  increase  of  temperature,  P  is 
expected  to  show  a  negative  temperature  dependence.  Figure  8  actually 
shows  P^  for  those  rotational  transitions  for  which  Aj^Q^  =  Aj^-Q  -  +  1. 

The  transitions  Aj^-Q  =  -  1  and  Aj^^  =  +  1  make  a  contribution  to  P^  that 
is  20  times  smaller.  In  the  second  order  calculation,  the  virtual  state  may 
be  an  excited  electronic  state  of  one  molecule  (Debye  induction  interaction) 
and/or  of  both  molecules  (dispersion  interaction).  In  Ref.  7,  it  was  shown 
that  the  contribution  to  the  rate  constan  uy  these  interactions  is  about  two 
orders  of  magnitude  smaller  than  the  reaction  paths  in  which  virtual  states 
are  excited  vibrational  levels  of  the  ground  electronic  states.  In  addition, 
following  Ref.  7,  we  will  neglect  those  virtual  states  in  which  both  the 
molecules  are  either  vibrationally  excited  or  have  no  vibrational  excitation. 
Figure  9  shows  three  sets  of  Feynman  diagrams  which  have  to  be  added  to 
obtain  the  cross-section  for  the  energy  transfer.  The  final  rotational  slates 
in  Fig.  9a  differ  from  those  in  Figs.  9b  and  9c  and  therefore  the  cross- 
sections  obtained  from  9a  add  to  those  obtained  from  diagrams  9b  and  9c. 

On  the  other  hand,  9b  and  9c  can  connect  the  same  initial  and  final  states 
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-3 


T  *K 


1  >c.  H  <■  0>  <  (>  |)ip.i'<-.<!ipoU>  transition  probability  vs.  temperature. 

I'  "  r  •>>>'.  v'O,  an<|  O  A  i  ■  +  1. 
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HS3S  1 


Fig.  9  CO^/CO:  Feynman  diagrams.  The  multipole  moments  connecting 
initial  states  with  a  virtual  state  are  shown  on  the  left  vertex, 
multipole  moments  connecting  final  states  with  a  virtual  state  are 
shown  on  the  right  vertex  (Q  =  quadrapole,  d  =  dipole). 
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B5709 


■  ■  •  — •••  |  •  ■  •  •  •  J . .  i  »»*  t  i»|/v  •>  »  Miwmv  vu  v>  v/  (  i  »  v  i  (  v  v  <  y  | 

itiul  priim's  denote  tln-ir  derivatives  with  respect  to  intermit- lea  r  distance 

evaluated  at  the  equi  I  il>  rinm  internuelear  distanee.  The  calculations  ***’  ** 

show  that  (J  and  (J*  have  opposite  signs.  The  calculations  of  CO  dipole 

moment  cited  earlier  give  the  sign  (CO  )  opposite  to  that  observed 
Z  1 

experimentally.  I'his  is  because  ot  small  dipole  moment  of  CO.  Minimi¬ 
zation  ot  nicrjiy  gives  a  slightly  distorted  charge  distribution  yielding 
opposite  signs.  If  we  accept  the  result  arrived  at  by  calculation***  that  at 
large  distance  the  sign  of  tin*  dipole  moment  is  C^O  ,  one  is  led  to  the 
c  om  lusion  that  dipole  moment  anti  its  derivative  have  opposite  signs.  In 
this  case  then,  the  interference  between  diagrams  9b  anil  9c  is  constructive, 
[n  call  ulating  the  contribution  of  diagram  9c  to  the  rate  constant  we  are 
again  confronted  with  an  uncertain  value  of  O'.  Fortunately,  the  fact  that 
CO  has  a  small  value  of  dipole  moment  is  of  help  here.  The  relative 
lontribution  of  l  ie.  ■»!>  and  ‘»c  to  the  scattering  amplitude  is 

*  1 

F  ir  the  scattering  amplitude  due  to  diagram  9c  to  be  about  0.  1  that  of  9b 

)  7 

1  Q  0  •  Z.  s  x  Id  “  .  I’his  number  is  approximately  Z.  5  times  the 
corresponding  matrix  element  for  N,  and  about  4  times  larger  than  value 
obtained  using  molecular  structure  i  ulculations.  '  *  Thus  the  inclusion 

i>:  diagrams  \  is  expet  ted  to  give  e  ros  s  -  sections  which  are  at  most  about 
1  i  larger  than  tin-  one  obtained  by  neglecting  it.  Since  other  assumptions 
n  the  ,  a  1 1  u  hit  ion  a  re  greater  sou  r  i  e  ni  error,  wt*  will  neglect  diagram  9  c . 


0  |  n  1  •  0  Q  0  • 

o  [L 1  <>'-x  rzTo  • 


0  I  Q I  0  ■  1  I  <  0  I  a  1 

*r" ,-r  - — t*~,  since  4 — — 


i  x  1  Q  0 


<0!pl0M 
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P  _•  5/T2  d'14  C2(246 


n,,2(f  H,^s6„.8,xz1 

VJ2  ' 


(35) 


(j,  .  j7.  *E)  • 


,u>\2 


where  (Q'  )  is  the  absolute  square  of  the  dipole  moment  matrix  element 

between  the  two  vibrational  levels,  and  is  the  expectation  value  of  the 
quadrupole  moment  of  molecule  { i ) .  The  expectation  values  of  CO^  and  CO 
in  the  excited  vibrational  state  are  taken  equal  to  their  values  in  the  ground 
vibrational  level. 


(36) 


2  +14 

(x,  8)  =  v  d 
b 

'+6 

z  (- 

m  =  6 

•  00 

where  x  =  ■  -  u>d 


(•*) 


1/2 


,  where  u>  =  2 n  c~  is  the  absolute  value 


of  the  energy  transferred  from  internal  degrees  of  freedom  to  translation. 
7,j  depends  upon  the  initial  and  final  rotational  levels  and  upon  the  energy 
mismatch  in  the  following  way: 

(i)  jj  =  +  3,  i  =  1, 2 

/  l.5(ji+l)  (j.  +  2)  (J.  +  3)  \fj-y.  -1 

Z,  <VJ2'  AE*  °\l=l,2  (2J.+  1)  7ZJTTJT  (2j.  +  5)  y[Lj"1>2  2Bi  <2Ji+3)]  - 


[ 


AE  - 


2  B 


i=  1,  2 


-In  2 

i(V»]  ] 


(37) 
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(u)  ->i  'l  '  1  •  '2  'j  '  5 


/ 


I.i,  (i,  ‘  n  (jj  U)  l.M.i,  M)  (j,  1 2)  (j,  f  3) 


'2  w  '-'2 


i  r?j  ,""'n 


ITJjr  n  (2j2  4  tt  (2jjnrr 


(  S  il':  Ujj  »  i) 

'  i  1,2 


7ZJ 


7TTT  ^TJ“T5-; 


(^K-  ‘  «■  (^fri  )  (AE  +  2B|j|-2B2<j2M))’1)] 


v  =  I  IF  -  2Bj  j  L  (j,  f  1)  -  6  Bzj2  (j2f  2)1 


(38) 


Figure  10  shows  a  graph  of  P  vs  T  for  cases  (i)  and  (ii).  Figure  11 
shows  the  P  due  to  first  order  and  second  order  process  and  their  sum. 
Also  shown  are  the  experimental  points  obtained  by  Fig.  6  and  Rapp's 
calculation  for  Fq.  (30). 
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(COLLISIONS) 


Fig.  10  COo/CO:  Various  second-order  contributions  to  transition 
probability. 

L  A  3co  =  Ajco2  =  +  3- 

2-  A  3co2  =  +  l'  ^CO  =  +  3* 

3-  ^-'co  ~  +  a,'co2  =  +  3- 

4.  Sum  of  1,2,  and  3. 
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long  range  forces  are  within  a  factor  of  two  off  from  the  experimental 
data.  In  addition,  both  show  approximately  the  same  temperature 
dependence  as  the  experimentally  observed  data.  The  long  range  forces 
calculation  does  not  contain  any  adjustable  parameters  and  its  validity  has 
been  established  by  studying  energy  transfer  from  N^O  to  ('.  ( Kq .  (31)  ). 

The  short  range  force  calculation  used  here  is  not  directly  applicable  to 
CO^-CO  energy  transfer  process  but  was  taken  from  a  calculation  of  V-V 
energy  transfer  between  two  diatomic  molecules  with  similar  energy 
mismatch.  Since  short  range  calculations  are  almost  always  adjusted  for 
the  range  parameter,  it  was  not  considered  advisable  to  invest  time  in 
another  short  range  force  calculation.  It  is,  however,  quite  concc  ./able 
that  Rapp's  calculation  overestimates  the  short  range  contribution  bv  a 
factor  of  4.  This,  then,  added  to  long  range  force  cross-section  would 
give  a  reasonable  fit  to  the  experimental  data.  This  explanation,  however, 
suffers  from  one  draw  back.  Why  should  the  contribution  to  the  cross- 

section  of  the  short  range  forces  to  CO^-CO  energy  transfer  process  at 

o  -  4  -5 

900  K  be  about  10  and  for  CO^-N^  reaction  less  about  10 

Another  possibility  is  that  short  range  forces  in  CO^-CO  energy 
transfer  process  do  not  play  a  role  for  the  temperatures  investigated  here. 
Short  range  forces  certainly  do  not  play  any  role  up  to  1000°K  for  CO^-N^ 
reaction.  So  this  possibility  is  not  obviously  incorrect.  The  difference 
between  the  long  range  calculation  and  experimentally  observed  data  may 
be  due  to  the  neglect  of  some  of  the  Feyman  diagrams,  for  example  those 
which  lead  to  A  =  Aj^  =  +  1.  This  possibility  is  currently  being  investigated. 
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